The source of the Black Death in fourteenth-century central Eurasia.

The origin of the medieval Black Death pandemic (AD 1346-1353) has been a topic of continuous investigation because of the pandemic's extensive demographic impact and long-lasting consequences1,2. Until now, the most debated archaeological evidence potentially associated with the pandemic's initiation derives from cemeteries located near Lake Issyk-Kul of modern-day Kyrgyzstan1,3-9. These sites are thought to have housed victims of a fourteenth-century epidemic as tombstone inscriptions directly dated to 1338-1339 state 'pestilence' as the cause of death for the buried individuals9. Here we report ancient DNA data from seven individuals exhumed from two of these cemeteries, Kara-Djigach and Burana. Our synthesis of archaeological, historical and ancient genomic data shows a clear involvement of the plague bacterium Yersinia pestis in this epidemic event. Two reconstructed ancient Y. pestis genomes represent a single strain and are identified as the most recent common ancestor of a major diversification commonly associated with the pandemic's emergence, here dated to the first half of the fourteenth century. Comparisons with present-day diversity from Y. pestis reservoirs in the extended Tian Shan region support a local emergence of the recovered ancient strain. Through multiple lines of evidence, our data support an early fourteenth-century source of the second plague pandemic in central Eurasia.

Stone Age Yersinia pestis genomes shed light on the early evolution, diversity, and ecology of plague.

The bacterial pathogen Yersinia pestis gave rise to devastating outbreaks throughout human history, and ancient DNA evidence has shown it afflicted human populations as far back as the Neolithic. Y. pestis genomes recovered from the Eurasian Late Neolithic/Early Bronze Age (LNBA) period have uncovered key evolutionary steps that led to its emergence from a Yersinia pseudotuberculosis-like progenitor; however, the number of reconstructed LNBA genomes are too few to explore its diversity during this critical period of development. Here, we present 17 Y. pestis genomes dating to 5,000 to 2,500 y BP from a wide geographic expanse across Eurasia. This increased dataset enabled us to explore correlations between temporal, geographical, and genetic distance. Our results suggest a nonflea-adapted and potentially extinct single lineage that persisted over millennia without significant parallel diversification, accompanied by rapid dispersal across continents throughout this period, a trend not observed in other pathogens for which ancient genomes are available. A stepwise pattern of gene loss provides further clues on its early evolution and potential adaptation. We also discover the presence of the flea-adapted form of Y. pestis in Bronze Age Iberia, previously only identified in in the Caucasus and the Volga regions, suggesting a much wider geographic spread of this form of Y. pestis. Together, these data reveal the dynamic nature of plague's formative years in terms of its early evolution and ecology.

Bacteriophages fEV-1 and fD1 Infect Yersinia pestis.

Bacteriophages vB_YpeM_fEV-1 (fEV-1) and vB_YpeM_fD1 (fD1) were isolated from incoming sewage water samples in Turku, Finland, using Yersinia pestis strains EV76 and KIM D27 as enrichment hosts, respectively. Genomic analysis and transmission electron microscopy established that fEV-1 is a novel type of dwarf myovirus, while fD1 is a T4-like myovirus. The genome sizes are 38 and 167 kb, respectively. To date, the morphology and genome sequences of some dwarf myoviruses have been described; however, a proteome characterization such as the one presented here, has currently been lacking for this group of viruses. Notably, fEV-1 is the first dwarf myovirus described for Y. pestis. The host range of fEV-1 was restricted strictly to Y. pestis strains, while that of fD1 also included other members of Enterobacterales such as Escherichia coli and Yersinia pseudotuberculosis. In this study, we present the life cycles, genomes, and proteomes of two Yersinia myoviruses, fEV-1 and fD1.

Yersinia pestis: the Natural History of Plague.

The Gram-negative bacterium Yersinia pestis is responsible for deadly plague, a zoonotic disease established in stable foci in the Americas, Africa, and Eurasia. Its persistence in the environment relies on the subtle balance between Y. pestis-contaminated soils, burrowing and nonburrowing mammals exhibiting variable degrees of plague susceptibility, and their associated fleas. Transmission from one host to another relies mainly on infected flea bites, inducing typical painful, enlarged lymph nodes referred to as buboes, followed by septicemic dissemination of the pathogen. In contrast, droplet inhalation after close contact with infected mammals induces primary pneumonic plague. Finally, the rarely reported consumption of contaminated raw meat causes pharyngeal and gastrointestinal plague. Point-of-care diagnosis, early antibiotic treatment, and confinement measures contribute to outbreak control despite residual mortality. Mandatory primary prevention relies on the active surveillance of established plague foci and ectoparasite control. Plague is acknowledged to have infected human populations for at least 5,000 years in Eurasia. Y. pestis genomes recovered from affected archaeological sites have suggested clonal evolution from a common ancestor shared with the closely related enteric pathogen Yersinia pseudotuberculosis and have indicated that ymt gene acquisition during the Bronze Age conferred Y. pestis with ectoparasite transmissibility while maintaining its enteric transmissibility. Three historic pandemics, starting in 541 AD and continuing until today, have been described. At present, the third pandemic has become largely quiescent, with hundreds of human cases being reported mainly in a few impoverished African countries, where zoonotic plague is mostly transmitted to people by rodent-associated flea bites.

Yersinia pestis Plasminogen Activator.

The Gram-negative bacterium Yersinia pestis causes plague, a fatal flea-borne anthropozoonosis, which can progress to aerosol-transmitted pneumonia. Y. pestis overcomes the innate immunity of its host thanks to many pathogenicity factors, including plasminogen activator, Pla. This factor is a broad-spectrum outer membrane protease also acting as adhesin and invasin. Y. pestis uses Pla adhesion and proteolytic capacity to manipulate the fibrinolytic cascade and immune system to produce bacteremia necessary for pathogen transmission via fleabite or aerosols. Because of microevolution, Y. pestis invasiveness has increased significantly after a single amino-acid substitution (I259T) in Pla of one of the oldest Y. pestis phylogenetic groups. This mutation caused a better ability to activate plasminogen. In paradox with its fibrinolytic activity, Pla cleaves and inactivates the tissue factor pathway inhibitor (TFPI), a key inhibitor of the coagulation cascade. This function in the plague remains enigmatic. Pla (or pla) had been used as a specific marker of Y. pestis, but its solitary detection is no longer valid as this gene is present in other species of Enterobacteriaceae. Though recovering hosts generate anti-Pla antibodies, Pla is not a good subunit vaccine. However, its deletion increases the safety of attenuated Y. pestis strains, providing a means to generate a safe live plague vaccine.

A genomic and historical synthesis of plague in 18th century Eurasia.

Plague continued to afflict Europe for more than five centuries after the Black Death. Yet, by the 17th century, the dynamics of plague had changed, leading to its slow decline in Western Europe over the subsequent 200 y, a period for which only one genome was previously available. Using a multidisciplinary approach, combining genomic and historical data, we assembled Y. pestis genomes from nine individuals covering four Eurasian sites and placed them into an historical context within the established phylogeny. CHE1 (Chechnya, Russia, 18th century) is now the latest Second Plague Pandemic genome and the first non-European sample in the post-Black Death lineage. Its placement in the phylogeny and our synthesis point toward the existence of an extra-European reservoir feeding plague into Western Europe in multiple waves. By considering socioeconomic, ecological, and climatic factors we highlight the importance of a noneurocentric approach for the discussion on Second Plague Pandemic dynamics in Europe.

New ancient Eastern European Yersinia pestis genomes illuminate the dispersal of plague in Europe.

Yersinia pestis, the causative agent of plague, has been prevalent among humans for at least 5000 years, being accountable for several devastating epidemics in history, including the Black Death. Analyses of the genetic diversity of ancient strains of Y. pestis have shed light on the mechanisms of evolution and the spread of plague in Europe. However, many questions regarding the origins of the pathogen and its long persistence in Europe are still unresolved, especially during the late medieval time period. To address this, we present four newly assembled Y. pestis genomes from Eastern Europe (Poland and Southern Russia), dating from the fifteenth to eighteenth century AD. The analysis of polymorphisms in these genomes and their phylogenetic relationships with other ancient and modern Y. pestis strains may suggest several independent introductions of plague into Eastern Europe or its persistence in different reservoirs. Furthermore, with the reconstruction of a partial Y. pestis genome from rat skeletal remains found in a Polish ossuary, we were able to identify a potential animal reservoir in late medieval Europe. Overall, our results add new information concerning Y. pestis transmission and its evolutionary history in Eastern Europe. This article is part of the theme issue 'Insights into health and disease from ancient biomolecules'.

Pentaplex real-time PCR for differential detection of Yersinia pestis and Y. pseudotuberculosis and application for testing fleas collected during plague epizootics.

Upon acquiring two unique plasmids (pMT1 and pPCP1) and genome rearrangement during the evolution from Yersinia pseudotuberculosis, the plague causative agent Y. pestis is closely related to Y. pseudotuberculosis genetically but became highly virulent. We developed a pentaplex real-time PCR assay that not only detects both Yersinia species but also differentiates Y. pestis strains regarding their plasmid profiles. The five targets used were Y. pestis-specific ypo2088, caf1, and pst located on the chromosome, plasmids pMT1 and pPCP1, respectively; Y. pseudotuberculosis-specific chromosomal gene opgG; and 18S ribosomal RNA gene as an internal control for flea DNA. All targets showed 100% specificity and high sensitivity with limits of detection ranging from 1 fg to 100 fg, with Y. pestis-specific pst as the most sensitive target. Using the assay, Y. pestis strains were differentiated 100% by their known plasmid profiles. Testing Y. pestis and Y. pseudotuberculosis-spiked flea DNA showed there is no interference from flea DNA on the amplification of targeted genes. Finally, we applied the assay for testing 102 fleas collected from prairie dog burrows where prairie dog die-off was reported months before flea collection. All flea DNA was amplified by 18S rRNA; no Y. pseudotuberculosis was detected; one flea was positive for all Y. pestis-specific targets, confirming local Y. pestis transmission. Our results indicated the assay is sensitive and specific for the detection and differentiation of Y. pestis and Y. pseudotuberculosis. The assay can be used in field investigations for the rapid identification of the plague causative agent.

The EnteroBase user's guide, with case studies on Salmonella transmissions, Yersinia pestis phylogeny, and Escherichia core genomic diversity.

EnteroBase is an integrated software environment that supports the identification of global population structures within several bacterial genera that include pathogens. Here, we provide an overview of how EnteroBase works, what it can do, and its future prospects. EnteroBase has currently assembled more than 300,000 genomes from Illumina short reads from Salmonella, Escherichia, Yersinia, Clostridioides, Helicobacter, Vibrio, and Moraxella and genotyped those assemblies by core genome multilocus sequence typing (cgMLST). Hierarchical clustering of cgMLST sequence types allows mapping a new bacterial strain to predefined population structures at multiple levels of resolution within a few hours after uploading its short reads. Case Study 1 illustrates this process for local transmissions of Salmonella enterica serovar Agama between neighboring social groups of badgers and humans. EnteroBase also supports single nucleotide polymorphism (SNP) calls from both genomic assemblies and after extraction from metagenomic sequences, as illustrated by Case Study 2 which summarizes the microevolution of Yersinia pestis over the last 5000 years of pandemic plague. EnteroBase can also provide a global overview of the genomic diversity within an entire genus, as illustrated by Case Study 3, which presents a novel, global overview of the population structure of all of the species, subspecies, and clades within Escherichia.

A single introduction of Yersinia pestis to Brazil during the 3rd plague pandemic.

Yersinia pestis was introduced to Brazil during the third plague pandemic and currently exists in several recognized foci. There is currently limited available phylogeographic data regarding Y. pestis in Brazil. We generated whole genome sequences for 411 Y. pestis strains from six Brazilian foci to investigate the phylogeography of Y. pestis in Brazil; these strains were isolated from 1966 to 1997. All 411 strains were assigned to a single monophyletic clade within the 1.ORI population, indicating a single Y. pestis introduction was responsible for the successful establishment of endemic foci in Brazil. There was a moderate level of genomic diversity but little population structure among the 411 Brazilian Y. pestis strains, consistent with a radial expansion wherein Y. pestis spread rapidly from the coast to the interior of Brazil and became ecologically established. Overall, there were no strong spatial or temporal patterns among the Brazilian strains. However, strains from the same focus tended to be more closely related and strains isolated from foci closer to the coast tended to fall in more basal positions in the whole genome phylogeny than strains from more interior foci. Overall, the patterns observed in Brazil are similar to other locations affected during the 3rd plague pandemic such as in North America and Madagascar.

Emergence and Spread of Basal Lineages of Yersinia pestis during the Neolithic Decline.

Between 5,000 and 6,000 years ago, many Neolithic societies declined throughout western Eurasia due to a combination of factors that are still largely debated. Here, we report the discovery and genome reconstruction of Yersinia pestis, the etiological agent of plague, in Neolithic farmers in Sweden, pre-dating and basal to all modern and ancient known strains of this pathogen. We investigated the history of this strain by combining phylogenetic and molecular clock analyses of the bacterial genome, detailed archaeological information, and genomic analyses from infected individuals and hundreds of ancient human samples across Eurasia. These analyses revealed that multiple and independent lineages of Y. pestis branched and expanded across Eurasia during the Neolithic decline, spreading most likely through early trade networks rather than massive human migrations. Our results are consistent with the existence of a prehistoric plague pandemic that likely contributed to the decay of Neolithic populations in Europe.

Integrative approach using Yersinia pestis genomes to revisit the historical landscape of plague during the Medieval Period.

Over the last few years, genomic studies on Yersinia pestis, the causative agent of all known plague epidemics, have considerably increased in numbers, spanning a period of about 5,000 y. Nonetheless, questions concerning historical reservoirs and routes of transmission remain open. Here, we present and describe five genomes from the second half of the 14th century and reconstruct the evolutionary history of Y. pestis by reanalyzing previously published genomes and by building a comprehensive phylogeny focused on strains attributed to the Second Plague Pandemic (14th to 18th century). Corroborated by historical and ecological evidence, the presented phylogeny, which includes our Y. pestis genomes, could support the hypothesis of an entry of plague into Western European ports through distinct waves of introduction during the Medieval Period, possibly by means of fur trade routes, as well as the recirculation of plague within the human population via trade routes and human movement.

[Source tracing of the Yersinia pestis strains isolated from Heqing county, Yunnan province in 2017].

Objective: To understand the genotype of the Yersinia (Y.) pestis strains isolated from Heqing county, Yunnan province in 2017 and provide evidence for the prevention and control of plague in this area. Methods: Ten Y. pestis strains isolated from Heqing were typed by the detections of different region (DFR) and clustered regularly interspaced short palindromic repeats (CRISPRs) as well as multiple-locus variable-number tandem repeat analysis (MLVA). And the results were compared with those of the 93 Y. pestis strains from the adjacent plague foci of Heqing obtained from the established database for clustering analysis. Results: The results showed that Heqing strains had the same type of DFR (Genomovar 05) and CRISPRs (Cluster Ca7, Type 22) with isolates from the plague focus in Lijiang. Heqing strains and Lijiang strains were in the same cluster in MST and only VNTR loci N2117 and M23 of Heqing strains were different from that of Lijiang strains. Conclusion: The Y. pestis strains isolated from Heqing in 2017 were highly homogenous with the strains isolated from wild rodents in plague focus in Lijiang, and Heqing plague might be the result of further southward spread of Lijiang plague.

Scaffolding of Ancient Contigs and Ancestral Reconstruction in a Phylogenetic Framework.

Ancestral genome reconstruction is an important task to analyze the evolution of genomes. Recent progress in sequencing ancient DNA led to the publication of so-called paleogenomes and allows the integration of this sequencing data in genome evolution analysis. However, the de novo assembly of ancient genomes is usually fragmented due to DNA degradation over time among others. Integrated phylogenetic assembly addresses the issue of genome fragmentation in the ancient DNA assembly while aiming to improve the reconstruction of all ancient genomes in the phylogeny simultaneously. The fragmented assembly of the ancient genome can be represented as an assembly graph, indicating contradicting ordering information of contigs. In this setting, our approach is to compare the ancient data with extant finished genomes. We generalize a reconstruction approach minimizing the Single-Cut-or-Join rearrangement distance towards multifurcating trees and include edge lengths to improve the reconstruction in practice. This results in a polynomial time algorithm that includes additional ancient DNA data at one node in the tree, resulting in consistent reconstructions of ancestral genomes.

Analysis of 3800-year-old Yersinia pestis genomes suggests Bronze Age origin for bubonic plague.

The origin of Yersinia pestis and the early stages of its evolution are fundamental subjects of investigation given its high virulence and mortality that resulted from past pandemics. Although the earliest evidence of Y. pestis infections in humans has been identified in Late Neolithic/Bronze Age Eurasia (LNBA 5000-3500y BP), these strains lack key genetic components required for flea adaptation, thus making their mode of transmission and disease presentation in humans unclear. Here, we reconstruct ancient Y. pestis genomes from individuals associated with the Late Bronze Age period (~3800 BP) in the Samara region of modern-day Russia. We show clear distinctions between our new strains and the LNBA lineage, and suggest that the full ability for flea-mediated transmission causing bubonic plague evolved more than 1000 years earlier than previously suggested. Finally, we propose that several Y. pestis lineages were established during the Bronze Age, some of which persist to the present day.

Ten years of surveillance of the Yulong plague focus in China and the molecular typing and source tracing of the isolates.

Plague, caused by Yersinia pestis, was classified as a reemerging infectious disease by the World Health Organization. The five human pneumonic plague cases in Yulong County in 2005 gave rise to the discovery of a Yulong plague focus in Yunnan province, China. Thereafter, continuous wild rodent plague (sylvatic plague) was identified as the main plague reservoir of this focus. In this study, the epizootics in Yulong focus were described, and three molecular typing methods, including the different region (DFR) analysis, clustered regularly interspaced short palindromic repeats (CRISPRs), and the multiple-locus variable number of tandem repeats (VNTR) analysis (MLVA) (14+12), were used for the molecular typing and source tracing of Y. pestis isolates in the Yulong plague focus. Simultaneously, several isolates from the vicinity of Yunnan were used as controls. The results showed that during the 10-year period from 2006 to 2016, an animal plague epidemic occurred in 6 of those years, and 5 villages underwent an animal plague epidemic within a 30-km2 area of the Yulong plague focus. Searching for dead mice was the most effective monitoring method in this plague focus. No positive sample has been found in 6937 captured live rodents thus far, suggesting that the virulence of strains in the Yulong plague focus is stronger and the survival time of mice is shorter after infection. Strains from Lijiang, Sichuan and Tibet were of the same complex based on a typing analysis of DFR and CRISPR. The genetic relationship of Y. pestis illustrated by MLVA "14+12" demonstrates that Tibet and Sichuan strains evolved from the strains 1.IN2 (Qinghai, 1970 and Tibet, 1976), and Lijiang strains are closer to Batang strains (Batang County in Sichuan province, 2011, Himalaya marmot plague foci) in terms of genetic or phylogenic relationships. In conclusion, we have a deeper understanding of this new plague focus throughout this study, which provides a basis for effective prevention and control.

Yersinia pestis detection by loop-mediated isothermal amplification combined with magnetic bead capture of DNA

ABSTRACT We developed a loop-mediated isothermal amplification (LAMP) assay for the detection of Y. pestis by targeting the 3a sequence on chromosome. All 11 species of the genus Yersinia were used to evaluate the specificity of LAMP and PCR, demonstrating that the primers had a high level of specificity. The sensitivity of LAMP or PCR was 2.3 or 23 CFU for pure culture, whereas 2.3 × 104 or 2.3 × 106 CFU for simulated spleen and lung samples. For simulated liver samples, the sensitivity of LAMP was 2.3 × 106 CFU, but PCR was negative at the level of 2.3 × 107 CFU. After simulated spleen and lung samples were treated with magnetic beads, the sensitivity of LAMP or PCR was 2.3 × 103 or 2.3 × 106 CFU, whereas 2.3 × 105 or 2.3 × 107 CFU for magnetic bead-treated liver samples. These results indicated that some components in the tissues could inhibit LAMP and PCR, and liver tissue samples had a stronger inhibition to LAMP and PCR than spleen and lung tissue samples. LAMP has a higher sensitivity than PCR, and magnetic bead capture of DNAs could remarkably increase the sensitivity of LAMP. LAMP is a simple, rapid and sensitive assay suitable for application in the field or poverty areas.

Molecular epidemiological investigations of plague in Eastern Province of Zambia.

BACKGROUND: Plague is a flea-borne zoonotic and invasive disease caused by a gram negative coccobacillus bacterium called Yersinia pestis. Plague has caused three devastating pandemics globally namely: the Justinian, Black Death and Oriental plague. The disease in the Eastern Province of Zambia has been reported in Nyimba and Sinda Districts in the past 15 years. The aim of this study was to investigate the molecular epidemiology of plague in the two affected districts. Polymerase Chain Reaction (PCR), targeting Plasminogen activator gene (pla gene) of Y. pestis, was performed on suspected human bubo aspirates (n = 7), rodents (n = 216), shrews (n = 27) and fleas (n = 1494). Of these, one positive sample from each source or host was subjected to sequencing followed by phylogenetic analysis. RESULTS: The plasminogen activator gene (pla gene) of Y. pestis was detected in 42.8% bubo aspirates, 6.9% rodents, 3.7% shrew and 0.8% fleas. The fleas were from pigs (n = 4), goats (n = 5) and rodents (n = 3). The sequencing and phylogenetic analysis suggested that the pla gene of Y. pestis in Nyimba and Sinda was similar and the isolates demonstrated a high degree of evolutionary relationship with Antiqua strains from the Republic of Congo and Kenya. CONCLUSION: It can be concluded that pla gene of Y. pestis was present in various hosts in the two districts and the strains circulating in each district were similar and resembles those in the Republic of Congo and Kenya.

Yersinia pestis detection by loop-mediated isothermal amplification combined with magnetic bead capture of DNA.

We developed a loop-mediated isothermal amplification (LAMP) assay for the detection of Y. pestis by targeting the 3a sequence on chromosome. All 11 species of the genus Yersinia were used to evaluate the specificity of LAMP and PCR, demonstrating that the primers had a high level of specificity. The sensitivity of LAMP or PCR was 2.3 or 23CFU for pure culture, whereas 2.3×104 or 2.3×106CFU for simulated spleen and lung samples. For simulated liver samples, the sensitivity of LAMP was 2.3×106CFU, but PCR was negative at the level of 2.3×107CFU. After simulated spleen and lung samples were treated with magnetic beads, the sensitivity of LAMP or PCR was 2.3×103 or 2.3×106CFU, whereas 2.3×105 or 2.3×107CFU for magnetic bead-treated liver samples. These results indicated that some components in the tissues could inhibit LAMP and PCR, and liver tissue samples had a stronger inhibition to LAMP and PCR than spleen and lung tissue samples. LAMP has a higher sensitivity than PCR, and magnetic bead capture of DNAs could remarkably increase the sensitivity of LAMP. LAMP is a simple, rapid and sensitive assay suitable for application in the field or poverty areas.